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Abstract 

We present a new prescription for the resummation of the divergent series of perturbative cor- 
rections, due to soft gluon emission, to hard processes near threshold in perturbative QCD 
(threshold resummation). This prescription is based on Borel resummation, and contrary to the 
commonly used minimal prescription, it does not introduce a dependence of resummed physical 
observables on the kinematically unaccessible x — > region of parton distributions. We com- 
pare results for resummed deep-inelastic scattering obtained using the Borel prescription and 
the minimal prescription and exploit the comparison to discuss the ambiguities related to the 
resummation procedure. 
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The resummation of logarithmically enhanced contributions to hard processes near thresh- 
old [1,2], such as deep-inelastic scattering and Drell-Yan production at large values of the Bjorken 
X variable (or its analogue in the case of Drell-Yan) , is characterized by the fact that the effective 
scale of the process is a soft scale related to the emission process. This means that for a process 
with hard scale the resummation of large logs of 1 — a; effectively replaces the perturbative 
coupling as{Q^) with q;s(Q^(1 — x)). In the space of the variable N which is conjugate to x 
upon Mellin transformation, where the resummation is more naturally performed, the effective 
couphng is as{Q'^/N) and the soft limit x — > 1 corresponds to — > oo. This result, which 
has been understood long ago [3] on the basis of an analysis of evolution equations in the soft 
limit, and more recently in terms of effective theories [4], is a simple consequence of the fact [5] 
that in the soft limit cross sections only depend on x through the soft scale Q^{1 — x), so this 
dependence can be renormalization-group improved using standard techniques. 

As the scale decreases, the strong couphng increases and eventually it blows up at the Landau 
pole, so when 

x = xl = 1-— (1) 

resummed results diverge, and physical observables can be determined only by specifying a 
prescription to treat this divergence. A simple option, already discussed in ref. [3], is to perform 
the resummation in x space and cut off the phase space integration so that the dangerous 
X > xl region is excluded. The option which is more commonly used, however, is to perform 
the resummation in A^ space, and reconstruct the result in x space by Mcllin inversion. In this 
case, if the Mellin inversion is performed order by order in perturbation theory, the series of 
resummed x-spacc contributions diverges, and the problem is turned into that of summing a 
divergent series [6]. 

A commonly used way of treating this divergent series is the minimal prescription (MP) [6] , 
which, as we shall discuss in more detail, is based on the observation that the Mellin inversion 
integral of the resummed A^ space result exists if performed along a suitable contour. Further- 
more, the divergent series obtained from the order-by- order Mellin inversion is an asymptotic 
expansion of this integral. The minimal prescription, however, has the shortcoming that upon 
convolution the partonic cross section does not vanish in the unphysical x > 1 region, which 
implies that physical observables pick up a power-suppressed contribution from the unaccessible 
X — >■ region of parton distributions. 

In ref. [7] some of us suggested instead that the divergent series could be summed using 
the Borel method, and showed how this can be done at the leading logarithmic level for the 
logarithmic derivative of the resummed partonic cross section. Here we show how to perform 
this Borel resummation at any desired logarithmic order for any physical observable (such as, say, 
the DIS or Drell-Yan cross section at the hadronic level): namely, we give here a general Borel 
resummation prescription (BP). The availability of several resummation prescriptions is per se 
useful as a way of estimating the uncertainty of the resummation procedure. More interestingly, 
we will show that the Borel prescription solves the aforementioned problem of the minimal 
prescription. Indeed, the BP leads to a resummed partonic cross section which has the form of 
an x-space plus distribution, such as found at finite perturbative order, and thus gives physical 
observables by a convolution with parton distributions in the standard way. This is achieved 
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through the inclusion of a higher twist term in the resummed result. 

We will first summarize the properties of the resummed result and in particular the divergence 
of the resummed perturbative expansion. We will then describe the Borel resummation of 
the resummed partonic cross section, and specifically discuss its dependence on the choice of 
higher twist terms included in it. Finally we will compare the Borel prescription to the minimal 
prescription, and in particular compare results for physical observables obtained using either 
method. 

In order to understand the origin of the divergence of resummed results, let us consider first 
as an example the computation of the resummed leading log expression of 

where a is the Mellin transform 

a (^,as{fi^),N^ = j dxx^'^a ^^,as{fi^),x^ (3) 

of an observable a as{fi^),x^ such as the Drell-Yan cross section or a deep- inelastic structure 

function, computed at the parton level. In the soft limit, 7(as(Q^), A^) is computed up to terms 
which do not grow as iV — > oo, and at the leading logarithmic level, it is a function of as(Q^) In 
only. Explicitly, 



^LL{as{Q'),N)=g, / —asiQ'/n) = -^\n 
Ji n (Jo 



/I'll _ ll-i I 

(4) 



1 + alnl 



where gi is a constant, a = 1 for deep-inelastic scattering and a = 2 for Drell-Yan, we have used 
the explicit leading log form of as{Q^), 



Oisi^jJi^) „ 33 — 2nj 



1 + /?oa,(/i2)ln^ 127r 

and we have defined 

a = aPoasiQ'^). (6) 

Clearly, 'yLii'^siQ'^), N) has a branch cut on the positive real axis of the complex plane, 
starting at the Landau pole of eq. ([5]), 

iVL = ei (7) 

But if ■jLLicisiQ'^), N) were the Mellin transform of some function PLL{as{Q'^),x), it would be 
regular above some abscissa of convergence Nc, i.e. for all Re(A^) > Nc. Hence, 'jLiicisiQ'^), N) 
is not the Mellin transform of anything. However, to any finite fixed perturbative order M the 
inverse Mellin transform of 7(as(Q^), A^) is given by 

p(*^)(a.(g^),x) = -f^X]LJL«^ 1 / dNx-^ln'^; N > 0, (8) 
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where all Mellin inversion integrals can be computed exactly [see the appendix, eq. (1651) ]. It is 
easy to see that the limit of P^^^^asiQ"^), x) as M — >• cxd diverges. Indeed, if the limit existed, 
then one could interchange the sum over k and the integral over A^, but the sum over k is then 
the Taylor expansion of 7LL(as(Q^), N) eq. ([2]), which has finite radius of convergence |A^| < Nl, 
whereas the integral extends to infinity. 

In ref. [7] we have computed the divergent series eq. ([H]) explicitly and summed it d la Borel. 
The approach of that reference however exploits the explicit form of 7ll(«s(Q^), eq. (jlj), 
and in particular the fact that the integrand in eq. (jlj) has a simple pole. We now present a 
generalization of that method which reduces to it in the case of 7ll(q;s(Q^), A^), but can be 
applied to any resummed quantity. 

We start with a Mellin-space resummed quantity E, function of as{Q^), In and possibly 
other kinematical variables such as the rapidity, which we will not indicate explicitly. This 
resummed quantity is related to the partonic cross section <^ eq- ([3]), or a quantity 

derived from it such as 7 eq. ([2]). Now, in the soft limit a eq. ([3]) can be expanded as [6] 



a 



1^' 



ctq exp 



In — fifi 



a In — )+ 5(2 ( a In 



N 



+ ... 



(9) 

where (Tq is the Born level result. It is thus convenient to expand the generic resummed quantity 
E as 



M 



E(a,(Q2),L) = lim V /^^(^.(g^))^^ 



k=l 



where we have defined 



In- 



(10) 



(11) 



L = aln — = apoOsiQ'^) ... ^, 

with a^lQ"^) not necessarily given by its leading order expression. In the case of the computation 
of the partonic cross section, S (^^(Q^), L) is explicitly given by 



a [ ^,a,(/i2) 



=ao[l + S(«,(Q2),L)]. 



(12) 



For a generic resummed observable the series eq. (fTO!) has finite radius of convergence | L | < 1 
dictated by the location of the Landau pole. Hence the term-by-term inverse Mellin of the series 
eq. (fTOl) is divergent. The divergent series can be determined explicitly [7] [see eq. (l65l) of the 
appendix] to compute the inverse Mellin transform of S {as{Q'^),L) eq. but with M kept 
finite, 

r-TV+ioo TAT 

^-^S(a,(g2),L), M finite . 



N—ioo 



2m 



■ X 



(13) 



We get 



{as{Q'),x) = 

M 



^R^Hx) 



X 



+ 0[{l-xf] 



n=0 



M 

E 

k=n 



+ 0[{l-xf] 



(14) 
(15) 



3 



where we have defined 

c^ = {k+l)hk+i; ^ = \n{l-x), (16) 
A*^")(z) is the n-th derivative of 

A(.) . fiy, (17) 

O [(1 — x)^] denotes terms which are nonsingular in the hmit x 1, and for brevity we have 
omitted the exphcit dependence of the coefficients hk and of R^'' on asiQ"^)- The divergent series 
which we wish to sum is then 

R{x) = hm R^'^ix). (18) 

The divergence of R{x) can be removed by performing a Borel transform with respect to a, 
which gives 

71=0 k=n ^ ' 

The inner series has an infinite radius of convergence because its coefficients are factorially 
smaller than those of the series eq. (fTOj) . Because A(z) eq. (fT7|) is an entire function of 2, it is 
easy to show that this implies that the outer series is also convergent. Indeed, because the series 
^k^k^^ is convergent with the same radius as the series eq. (fTOl) . as /c — *• oo the coefficients Cfc 
are bounded by some constant K > 0, |cfe| < K. But this implies 

i^(^, x)\ < y -±—\wH'--\ = ifeM y \^—^\wr (20) 

' ^ ~ ^ n! ^(A;-n)!' ' ^ n\ ^ ^ 

n=0 k=n ^ ' n=0 

which converges because of the absolute convergence of the power series for l^{z) eq. (fTTI) . 
The original series can be recovered by inverting the Borel transform, 

/■oo 

R{x)= / dwe-'^R{w,x), (21) 

but the integral over w in eq. (1211) diverges at infinity: indeed, if we integrate the series term by 
term we recover the original divergent series eq. f|T5|) . We can cut off the singularity by extending 
the integral only up to some upper bound C . Because the series eq. (|T9|) converges uniformly in 
the interval < w < C, we can integrate term by term, with the result 

Rb[x,C) = / dwe-^R{w,x) (22) 



n— n h — n ^ ' 



(23) 



n=0 k=n 

where 



n=0 
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is the truncated gamma function. The series eq. fl23l) for Rb {x, C) has infinite radius of conver- 
gence, hke that for R{w,x) eq. f|T9|) . 

The Borel resummation of S {as{Q^), x) is obtained substituting the expression for Rb {x, C) 
eq. ( !23|) in eq. ( fMl) . Equation ( l23l) is not vey useful because it requires the evaluation of a 
double series. However, we will now show that the series can be summed through an integral 
representation which is not harder to evaluate numerically than the minimal prescription. Before 
doing this, let us discuss the properties of the Borel resummation. 

First, it is easy to see that the divergent series we started from R{x) eq. (fTSll is an asymptotic 
expansion of its Borel resummation Rb {x, C) eq. (!23l) . To this purpose, we note that R{x) and 
Rb {x, C) are related by 

RB{x,C) = Rix)-Rht{x,C), (25) 

where, using eq. ([2 



C) = e-« f A<") (1) f ( M c..-' t (I) 

n=0 k=n ^ ^ m=0 ' ^ ^ 



(26) 



Hence, Rm ~ e~^, so it vanishes faster than any power of as(Q^) as as(Q^) 0. It follows that 
the difference between Rb {x, C) and the sum of the first terms of R (x) is of order as{Q'^)^^^ ■, 
which proves that R {x) is an asymptotic expansion of Rb {x, C). 
Furthermore, note that using the expression for asiQ"^) 

asm = -7^[^ + 0{as{Q'))\ (27) 

we get 

e-^=(^^j [l + 0(a,(Q2))]. (28) 

This shows that cutting off the Borel inversion integral eq. (|2T1) at w = C is equivalent to 
including a twist-t contribution Rht{x,C), with 

2C , , 

t = 2 + — . 29 
a 

The divergence of the higher twist term then cancels that of the divergent series, leading to a 
finite result. The value of C should be chosen in such a way that no new, spurious higher twist 
terms are induced in physical observables. The choice C = a is minimal in that it corresponds 
to the inclusion of a twist-four term, i.e. a term of the first subleading twist. 

Let us now turn the Borel resummed expression Rb {x, C) eq. ( l23i) into a more useful form. 
In order to perform the sum over n, we introduce the Fourier transform 



r ^e-<''A(l + r/)e(r/ + r/o 



which satisfies 



(C) = / -^e-^-'A (1 + r/) e(r/ + r/o) (30) 



/oo 
dCe'^^A,, (C) = A(l + r^)0(r^ + r/o), (31) 
-oo 
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Figure 1: The Hankel contour used in eq. ([3; 



where Q{ri) is the Heaviside step function, and it is necessary to introduce a cutoff at 770 because 
the Fourier transform of the function /S.{z) does not exist. Rewriting Rb (x, C) eq. fl22fl9l) with 



we can perform the sum over n exphcitly: 



(32) 



Rb 



c 



dw e 



/oo °° 
fc=o 



(33) 



where it is sufficient to choose ?7o > to ensure that the result is independent of the choice of 

Now, we observe that the sum over k can be performed exphcitly if the factor of k\ in the 
denominator is removed. We do this as follows. First, we write l/k\ = l/r(fc + 1) and we use 
the Hankel representation of the Gamma function 



1 



Viz) 



H 



dt 

27ri' 



[-tr 



(34) 



where H is the Hankel contour shown in fig. [T] Furthermore, because the integrand in eq. ([3] 
doesn't have any other singularities in the complex plane besides the cut along the positive real 
axis, the integral along the Hankel contour H is equal to the integral along the contour Hi 
defined by 

ZH, = Re'', (35) 

with R 00 and e < 9 < In — e with e — > 0. If we substitute in eq. ( !33|) the expression eq. fl34|) 
with z = k + 1 and the integral over t performed along Hi we can integrate term by term over t 
the sum over k, because the contour Hi is always within the radius of convergence of the series 
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if R eq. is large enough. We get 



Rb 



dw e 



dC\o (0 



c 



— oo 
oo 



dwe ^ rfCA^,, (C) 



^-oo J Hi ^''^'^ 



dt 



k=0 



to 



where we have defined 



(36) 



(37) 



in terms of the function E (^^(Q^), L) eq. ( ITOl) . 

We can now remove the dependence of S' on ( through the change of variables 



w{l + <) 



(3^ 



whereby the contour Hi eq. (l35l) is mapped onto a contour f^i, which can be deformed back to 
the contour Hi for the new variable ^. The integral over ( can then be performed using eq. ( |3T1) 
with the result 



Rb (x, C) 



dw e 



Hi 



(l-x)"«A(l+0 S'(l/e). 



(39) 



The result eq. (139|1 can be already used as a resummation prescription. However, it may be 
more convenient to rewrite it directly in terms of the resummed observable S rather than its 
partial derivative. This is accomplished integrating by parts: 



Jo w Jh^ 



d^ d 
2^idi 



w 



ee-^«A(l + <) 



s(i/e) 



(40) 



where the surface term vanishes provided only the radius R eq. fl35l) of the contour in the ^ plane 
is large enough, because S(l/^) has a discontinuity along the negative real ^ axis that only 
extends from the origin up to the location of the Landau pole at = — 1. With straightforward 
manipulations we can rewrite eq. (140|) as 



Rb= [ ^\wiC,i,0 + - rdwW{w,i,0 

J Hi 2772 L a Jo 



s (1/0 



(41) 



where we have defined 

W{w,e,0 = we-^^^-^^''^A{l + wO =we-^(l-x)"'«A(l+0. (42) 
The Borel prescription for the resummation of the divergent series eq. (1151) consists of taking 

'RB{x,cy 



1 — X 



(43) 



with Rb{x, C) given by either of the equivalent expressions eq. (15^ or eq. (HTj) . The integrand of 
the ^ integral has a cut along the negative real ^ axis for —1 < ^ < 0, and it is regular elsewhere; 
the closed contour Hi encircles this cut. The value of C is related by eq. fl29|) to the twist of 
the contribution which is included in order to get a finite resummed result; the minimal choice 
is C = a, corresponding to the inclusion of a twist four term. 

Let us now briefly discuss some properties of the Borel resummed result eq. (H3l) and then 
compare it to the result obtained using the minimal prescription. First, let us determine it 
explicitly in the simplest case in which we take as resummed observable 

E(«,(Q2),l) =7LL(«s(g'),iV), (44) 
with 'jLL{cis{Q'^), N) given by eq. (jlj). In this case, it is convenient to use eq. (1371) . since 

s'(i/0 = -f TT7 (45) 

so the ^ integral is straightforward: 

RB{asiQ^),x,C) = dwe ^^Ail-w) 

Po Jo 



' A{l-w) (46) 



Jo 



A2 



Q2(l-x) 



which coincides with the result of ref. [7]. The next-to- leading log result can be analogously 
determined in closed form. 

The resummed result eq. ( l43l) has the form of a plus distribution, whose action on any test 
function /(x) leads to a finite result provided the numerator Rb{x, C) is integrable as a function 
of X between and 1. However, the explicit result eq. (H6!) suggests that this is the case only if 
C is not too large. Indeed, the integrand of eq. (!39|) is integrable over x as x ^ 1 only if 

Re«) > -1. (47) 

The path Hi must intersect the negative real axis at some ^ = < ~1 because of the cut up to 
the Landau pole at ^ = —1. Hence the condition becomes Re(w) < 1 which is violated whenever 
C>1. 

Nevertheless, for all C the action of S (^^(Q^),^) is well defined by analytic continuation, 
and this allows its numerical implementation. Indeed, consider the action upon integration of 
E(as(Q^),x) on the test function 

t{x) = (1 - x)". (48) 

We get 



fdxJ:[as{Q^),x] r(x) = -i. r^^e-t E'(l/e) A(l + O^-. 
Jo 27r« Jo J i n + wi 

The integrand is regular at w = —n/C,, because, for any negative integer —n, 



(49) 



[z + n) 



(-l)"-i(n -!)![! + 0(2 + n)], (n > integer) (50) 
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hence the integral eq. fl49p exists for all C. This immediately implies that is 
a distribution which gives finite results when integrated over any test function r(x) which is 
analytic in the neighbourhood of x = 1. 

In practice, for numerical computations one can proceed as follows: the quantity of interest 
is typically a convolution of S {as{Q'^)^x) with a parton density q{x), of the form 



Rb{z,C) 



1 fx 



z 



q{x) 



q{x) 



(51) 



The first term on the right-hand side of this equation only leads to a convergent integral if 
Rb{z,C) is integrable, which in turn requires C < 1 as discussed above. However, we can 
rewrite this integral defining 



g{x,z) 



1-z 



1 /X 

-9 - 



z 



q{x) 



(52) 



as follows 



dz 



Rb{z,C) 



~^{~)~^^^'> dzRB{z,C)[g{x,z) - g{x,l)]+ g{x,l) j dzRBiz,C). 



(53) 



The second integral on the right-hand side of eq. fl53l) can be computed analytically using eq. fH9|) 
with n = 0, while the first integral is now convergent even if Rb{z, C) is not integrable, provided 
only Rb{z,C) ~ (1 — z)^ with b > —2, which now only requires C < 2. If Rb{z,C) is even 

Z->1 

more divergent around z = 1 one simply iterates the procedure. This allows one to choose an 
arbitrarily large value of C. 

We finally compare the Borel prescription eq. fH5]l to the minimal prescription. Consider first 
what happens when we apply either of them to a quantity whose inverse Mellin transform does 
exist, such as Tj{as{Q'^), L) eq. (fTOj) when M is kept finite. In such case, the minimal prescription 
simply gives this inverse Mellin transform. Taking for example eq. (11 01) with M = 1, hi = 1, 
i.e., T.^asi^Q'^), L) = L, the minimal prescription gives [see appendix, eq. (^5^] 



^l,MP, 



,x 



a 



Ini 



+ S{1 



x] 



(54) 



If instead we apply the Borel prescription, we get a result that differs from the inverse Mellin 
first, because terms which are either finite or zero as x — ^ 1 are neglected, and furthermore, 
because the higher twist correction eq. ( !26l) is included. In the previous example, this gives, 
instead of eq. ( l54l) . 

1 



.1,BP 



{osiQ ),x) = a 



1 — a; 



1-e 



(55) 



If one applies the MP to a function S(a;s(Q^), L) whose Mellin transform does not exist 
because of a branch cut from eq. ([7j), such as a typical resummed quantity, the ensuing 
x-space result T,^^{as{Q'^),x) does not vanish in the unphysical region a; > 1. It follows that a 
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physical observable a^^{x), computed combining a partonic cross section a^^{N) eq. with 
a parton distribution q{N) [with inverse Mellin q{x)], has the form 



i.e. it receives a contribution from the unphysical < y < x region of parton densities (see 
appendix B of ref. [6]), though it has been shown in ref. [6] that this contribution is power 
suppressed. Furthermore, there are practical difficulties in the construction of the x-space result 
iy^^{as{Q'^),x) which is needed e.g. if one wants to use a resummed result with x space parton 
distributions, related to the fact that the MP result for T,^^{as{Q'^),x) displays an oscillatory 
behaviour [6]. 

The main advantage of the Borel prescription result eq. (I43l) is that it gives directly T,^^ {as^Q"^), x) 
in X space, in the form of a plus distribution as those found order by order in perturbation theory. 
Physical observables are obtained from it by standard convolution with parton distributions in 
the physical region: 



This is accomplished by including power suppressed terms order by order in the physical region, 
as explicitly shown in eq. fl55|) . As already mentioned, it is convenient to choose C in such a way 
that these power suppressed terms combine with those which already appear at higher orders 
in the Wilson expansion. In fact, it has been argued in ref. [8] that in the large x limit the 
dominant higher twist contributions are those which mix upon renormalization with the leading 
twist. Be that as it may, with the minimal choice C = a eq. fl2^ the ambiguity introduced by 
the BP may be cancelled by an equal and opposite ambiguity from a conventional higher twist 
term, as already discussed in ref. [7]. 

A further advantage of the BP is that the non-logarithmically enhanced terms which are 
generated by the exact Mellin inversion of S(as((5^),L) can be included or excluded at will. 
Indeed, the computation of the exact Mellin inverse, as done in the MP, is not necessarily 
advantageous if the resummed S(q;s(Q^),L) is only computed in the large N limit to begin 
with. For instance, in the simple example considered above, the series of terms generated by 
the expansion of 1/ In ^ = 1/(1 — x) — 1/2 — (1 — x)/12 + . . . in eq. fl5^ does not necessarily 
provide a better approximation to the exact 0{as) expression of S than the purely logarithmic 
contribution 1/(1 — x) included in the BP result eq. ([S5D- This is to be contrasted with the case 
of specific classes of non-enhanced [9] or even suppressed [10] terms whose resummation might 
be advantageous. Now, in the BP it is possible to choose whether to perform the Mellin inversion 
exactly or in the large x limit, unlike in the MP where the inversion is always performed exactly. 
Indeed, it is easy to modify the BP in such a way that when applied to T,{as{Q^), L) eq. (fTOj) 
with finite M it coincides with its exact inverse Mellin up to higher twist terms. For this, it is 
sufficient to use the method described in the appendix to determine the Mellin inversion eq. fl65p 
exactly. In practice, it is sufficient to replace everywhere (1 — x) with In - in eqs. f flUfTSl) and in 
the final results eq. fl39|) or eq. fH2l) . 

The ambiguity in the resummation procedure can be estimated by comparing results obtained 
using the Borel and minimal prescriptions. In order for the comparison to be significant, we must 




(56) 




(57) 
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compare the convolution of the result with a test parton distribution: indeed, the resummed 
partonic quantity Il{as{Q'^),x) is a distribution, rather than a function proper. Furthermore, the 
different treatment of non logarithmically enhanced terms between BP and MP is only allowed 
in the region where the resummed logs are large: indeed, away from that region any resummed 
result must reduce to the fixed order. Hence, we must compare results matched to the fixed 
order. 

To this purpose, we have determined a matched result for a physical observable using the MP 
and BP. We consider the quark coefficient function for the deep-inelastic structure function F2 
in the MS scheme, (Q^Z/i^, N). We then determine the resummed expression for this 

quantity with //^ = Q^, up to the next-to-leading log level, i.e. we use eq. for a (1, as{Q^), N), 
including the contributions gi and g2, as given e.g. in ref. [11]. This is our resummed observable 
E(a,(Q2)^ L), to be matched to the standard 0(aj result for C|f (1, a^Q"^), N) [12]. We further 
take a model quark distribution given by 

q{x) = x-'/\l-xf; g(Ar) = (58) 

The minimal prescription is then constructed by computing 

F^^{Q',x) = I ^x-^ [Cff (l,a(Q^),iV) +E(a.(Q^),L) -S^^O(«.(g2),L)] g(iV), 

(59) 

where Cmp is the standard minimal prescription contour [6], T,{as{Q'^), L) is the resummed 
coefficient discussed above, and T}^^'^{as{Q'^),L) is its expansion up to order ^^(Q^), namely 

S(a,(g2), L) - S^LO(a,(Q2), L) = O {a\Q')) . (60) 

The Borel prescription is constructed computing 

F^^,{Q\x) = [A^(l,a(Q^),i/) + S(a.(Q2),y) -S^LO(a.(Q2),y)] q(-\ , (61) 

where C|f (1, a(Q2), x) is the inverse Mellin transform of {l.a^Q"^), N), ^:{a,{Q'^),x) is 
constructed from 'L{as{Q'^), L) using eq. (H3!l with Rb{x,C) given by eq. fHT!) and (for DIS) 
a = 1. We take C = 1, which corresponds to the inclusion of a twist-four term; the convolution 
integral in eq. (!6T!) can then be computed with one subtraction eq. (!53l) . Finally, T,^^^{as{Q'^),x) 
is the expansion of T,{as{Q'^),x) up to order ^^(Q^), 

f:{a,m,x) - S^LO(a,(g2), x) = O {a\Q')) . (62) 

The results obtained using the MP and the BP are compared to each other and to the fixed 
0(q;s(Q^)) result in Fig. [2l The structure function eqs. fl59|) . fl6Tl) is plotted as a function of x, 
normalized to the parton distribution eq. (!58l) : namely, we plot -^"^^t^- We take as(Q^) = 0.2. 
Note that q{x) vanishes very rapidly as x ^ 1. The comparison shows that the effect of the 
resummation is sizable for x ^ 0.6 and becomes of order 100% when x ^ 0.8, where however 
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Figure 2: Resummation of the quark coefficient function {Q"^ / fj? , a^fx^) , N) for the deep- 
inelastic structure function F2. The resummation is performed up to the next-to-leading log- 
arithmic level matched to the 0(q;s(Q^)) fixed order result. We plot as a function of x the 
structure function normalized to the parton distribution. The three curves are, top to bottom, 
the minimal prescription eq. fl59l) . the Borel prescription eq. fl6T|) and the fixed 0{as{Q'^)) result. 
We take ^.(Q^) = 0.2. 



^2(0;) is very small. Interestingly, while the percentage difference between the MP and BP 
tends to zero both as x — > 1 and x — > 0, in the intermediate region 0.6 ^ x ^ 0.8 where the 
resummation is important the two prescriptions lead to rather different results. 

We have checked that the replacement of (1 — x) — In ^ in the BP has a negligible effect. 
This ensures that the difference between the MP and BP is not due to a different treatment of 
non-logarithmically enhanced terms. Furthermore, we have verified that increasing the value of 
C eq. ([22D from one to 1.8 also has essentially no effect. This agrees with expectations based 
on the results of ref. [7]: there, it was found that Rb{x, C) is stable upon variations of C unless 
a; > eq. ([I]), so the same should hold for physical observables where the region of very large 
X shouldn't weigh too much. Finally, in Fig. [3] we repeat the same calculation but adding an 
extra logarithmic order in the resummed result and matching to the fixed order result computed 
at 0{a1). The difference between the fixed-order and resummed results is now smaller, as it 
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Figure 3: Same as figure 2, but with the resummation performed up to the next-to- next-to- 
leading logarithmic level and matched to the O(a^) fixed order result. 



ought to be, but the difference between MP and BP has not decreased, thereby showing that 
this difference is not compensated by the inclusion of higher logarithmic orders. 

We must conclude that the difference between the MP and the BP indicates that the ambi- 
guity in the resummation procedure is sizable: a fact which is rather well known in the context 
of transverse momentum resummation (see e.g. ref. [13]), but not equally obvious for threshold 
resummation. 

In summary, we have presented a new prescription for the resummation of the divergent series 
of logarithmically enhanced terms which is obtained from threshold resummation. The divergent 
series is summed through the Borel method, and the divergence in the Borel inversion integral 
is removed through the inclusion of a suitable higher twist term. This term can be chosen to 
be of any twist, but the minimal choice is to take it as a twist four contribution. We have 
described the practical implementation of this prescription and demonstrated its application to 
the threshold resummation of a deep-inelastic coefficient function, which we have compared to 
the commonly used minimal prescription. 

The Borel prescription and minimal prescription have somewhat complementary advantages 
and disadvantages: the minimal prescription is naturally implemented in N space, so it is easy 
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to use with A^-depcndcnt parton distributions. However, in x space the MP leads to partonic 
cross sections which do not vanish in the unphysical x > 1 region and its implementation is less 
straightforward. The Borel prescription directly gives an x space result which has the form of a 
plus distribution such as found in fixed order perturbative computations. However, its N space 
form can only be obtained by performing the Mellin transform numerically, and its convolution 
with a parton distribution must be determined by numerical integration. The Borel prescription 
also has the advantage that it is possible to control the inclusion of non-logarithmically enhanced 
terms in the resummation, but it has the disadvantage that it requires the inclusion of higher 
twist contributions. 

Comparison of results obtained using the Borel prescription and the minimal prescription 
suggests that the ambiguity in threshold resummation is sizable. The extension of this method 
to the case of resummation of transverse momentum distributions will be presented elsewhere. 

Acknowledgement: We thank Paolo Nason for illuminating discussions. S.F. acknowledges 
partial support from the Marie Curie Research Training Network HEPTOOLS under contract 
MRTN-CT-2006-035505. 



Appendix 

In ref. [5] we have determined the Mellin transform of any function of IniV to all orders in 
ln(l — x), up to terms which vanish as x — > 1 as a power of 1 — x. However, the exact Mellin 
transform can also be computed [7]. Indeed, the standard Euler integral representation of the 
Gamma function implies that 



dxx 



N-l 



In"-^ 



X 



Tin) {N-^ - 1) . 



(63) 



so 



-1 fN+ioo 

— _ dNx-^N-1 = A{r]) 



In-'-i - 



X 



+ d{l-x), 



+ 



where A (77) = It follows that the exact inverse Mellin transform of In*^ -h is 



N 



(64) 



1 pN+ioo 1 1 fN+ioo 

/ dNx-"" In^i- = TT- / dNx 
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hTi^ 
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In"-^ - 
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T]=0 

+ S{l-x) 



r]=0 



k 



n 



nA("-^^(l) ( Inln- 
V X 



k—n 



(65) 
+ 6koS{l-x), 



+ 



where in the last step we used A(")(0) = nA^" ^^1)- 
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